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O ■ Abstract 
O 

' We determine new constraints on the products of two i?-parity violating coupling 
constants from the NOMAD and CHORUS experiments on — > oscillations. 

CIhI We obtain improved results from (a) lepton flavour violating (LFV) meson decays, 

^ ' (b) LFV deep-inelastic scattering and (c) a combination of the two. 

> 

X 1 Introduction 



The superpotential of the MSSM+^p is obtained from the superpotential of the MSSM 
by adding the following terms (c./. ref.|l[|) 

+1 e^y^X'U,, U\^D\''d\'' + e^'n, L\H^,. (1) 

H, Q, L represent the left chiral S'f/(2)vt'-doublet superfields of the Higgses, the quarks 

and leptons; U,D,E represent the right chiral superfields of the w-type quarks, d-type 
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quarks and electron-type leptons respectively; a superscript denotes charge conju- 
gation; a,b and x,y,z are SU{2)w and SU{3)c indices; i,j,k and later also f,l,m are 
generation indices (summation over repeated indices is implied); 6^^ is the Kronecker sym- 
bol, symbolizes any tensor that is totally antisymmetric with respect to the exchange 
of any two indices, with e^^ " = 1. The coupling constants Xijk,X"ijk are antisymmetric 
with respect to the exchange of the first two/last two indices. The last term in eq.(|l|) 
can be rotated away utilizing a unitary field-redefinition [^. 

Good agreement between the SM theory and experiment gives stringent upper bounds 
on the extra 45 coupling constants Xijk, X'ijk and X"ijk, as well as on products thereof. 
For a list of references and the processes dealt with, see e.g. [0, ^, ^. The first systematic 
analysis of^p bounds was presented in 0, among other processes also utilizing neutrino 
deep-inelastic scattering. We look at this in considerably more detail here including also 
lepton flavour violating meson decays and including explicitly experimental effects. We 
also go beyond to allow products of couplings. 

Both processes can be tested by neutrino oscillation experiments. In this article we 
concentrate on the NOMAD and CHORUS [|] experiments at CERN, for several 
reasons. These two experiments looked for muon-neutrinos oscillating to tau-neutrinos, 
thus addressing the products of couplings for second and third generation leptons, which 
are less constrained than the first generation. The limits obtained by these experiments 
are the most stringent to date, thanks to the high statistics running. Lastly, although the 
two experiments are dismantled, the two collaborations are still active, so the numerical 
factors (see below) which are necessary to extract firm limits and which are not available 
in the published literature can still be evaluated. 

2 Lepton Flavour Violation at CHORUS & NOMAD 

CERN's SPS fires a proton beam on a Beryllium target. From the resulting jet the 
secondary mesons (mostly pions and kaons) are charge selected and left to decay into 
mainly muons and neutrinos. After the neutrinos have propagated some distance they 
are indirectly detected via neutrino deep-inelastic scattering (z/-DIS) on a nucleon in the 
target, producing a charged lepton plus other debris. The dominating SM process is 

p + Be =^ (tt or A ) — > fi + 

SM 
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resulting in an isolated muon and a jet of many hadrons (X). We are interested in 
detecting deviations from the SM via lepton flavour violating interactions, in particular 
those resulting in a final state tau. Allowing for at most one non-SM process, we have 
the following possibilities 

+ ^p. ..c 



(/) p + Be ^ (7r+ or ^ + 



SM 



v, + N ^ r + X, 



(//) p + Be ^ (7r+ or K+) ^ fi^ + ^, 



oscil. 



v, + N ^ r + X, 



[III) p + Be ^ {71+ or K+) ^ fi^ + i^^, 

SM 



u^ + N ^ T + X. 



Here we focus on the cases (/) and (///), as well as the combination of (/) and (///). 
We do not further discuss (//) since the case of neutrino oscillations has been treated 
extensively, see e.g. [0, |1T[. A rough estimate of the case (///) was made in |T^, which 



mainly focuses on other issues. The detector efficiencies, the spectrum of the neutrino 
beam and the convolution of the partonic cross section with the parton distribution 
functions (pdf's) were not taken into account. In addition, the coupling considered in 



12 1 also contributes to (/) and the combined effect must be taken into account, as we do 
in section ^ Our bound is a factor 3.3 better. 

We are interested in the theoretical prediction of the ratio of the number of tau- 
detections Afr over the number of muon-detections A/"^ , to be compared with rjexp , the 
experimental upper bound on the number of tau-detections per muon-detection. The 
latter is given by CHORUS and NOMAD from their neutrino oscillation analysis, of 
which we make explicit use. We assume that the new physics does not significantly affect 
the muon detection rate and obtain 

A/" ftp- decay _j_ j^J'oscil. _j_ j^J'flp-DIS ftp- decay _j_ j^J'oscil. _j_ j^J'ftp-DIS 

K ^ A/;f ^ + Af^''-'''"'^ + AT--'- + Af^^-""'' W 



< 



exp- 



(2) 



^ NOMAD and CHORUS are also sensitive to i?-parity violation via late decaying neutralinos |^ . 



Here is the number of observed leptons I which are due to the process Y . We shall 
here only analyze r's which originate either from a l^p meson decay or from ^p-DIS. 
The probability that a r is detected due to flavour violating Z/-DIS is given by 

Pa = X n 5 ^ • (3) 



E"1\.T / SM, all El, I i 

E E (cT +^^^^+^ + a, 

I=n,p 1=1 1=1 



SM, a// E„ 

N=n,p 1=1 j=l ^ " uJ 



Here = Pa{^'3jk^2ik^ ^kji^iks)- is the number of the nucleon-type N in the detector 
target and a denotes a cross section which has been integrated over the indicated parton 
distribution function and over the incoming neutrino energy distribution. The detector 
efficiency correction factors, x and Xqinqa^ty different for each physical process, 

Rp SM all 

XqinQaut SM ' ^ SM all E^ ' 

e denotes the efficiency of detecting a tau-neutrino with the experimental cuts used in 
the selection analysis, for each respective process. The ratio of efficiencies needs to 
be determined case by case for each considered analysis by generating a sample of fip 
interactions and applying the experimental cuts on it. We work with the target being 
isoscalar, which is approximately correct for both experiments [|^ |^. Thus the factors 
N^r in eq.@ can be taken outside the sums and cancel. 

In order to analyze the lepton flavour violating meson decays we define the fractions 

P (\'* \' \'* \ \ — n^fJ,C:U-r TT^fl^u^ , . 

^TrlAgifcAsife, Afcn^3fc2j = — i^; YSM ' 

^A'l^32fc^21fc> ^fcl2^3fc2j = '~7^r' ' ^ > 

The denominators are dominated by the SM decay. For each decay mode we have ex- 
plicitly indicated the possible relevant coupling constants. 

Furthermore, is the fraction of tau-neutrino events having a pion as a parent meson 
which survive the cuts applied for the oscillation analysis. The corresponding fraction 
originating from kaon decays is denoted and to a good approximation F^ = 1 — -F^. 
Using eqs.(^^ with the relevant coupling constants, and since the probabilities are 
small, we obtain from eq.(0) 

Fn ■ Pni^SlkKlki Kll^3k2) + ~ ^Sj ' PK{>^32kKlki Kl2^3k2) + 

Poscil. + -P(7(A3*fcA2jfc, A'^* jA2fc3) < Vexp, (7) 

with i,k = 1,2,3 and j = 1, 2. We assume that the three processes do not mutually 
cancel each other, and we can thus determine separate bounds. Furthermore, we assume 
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in turn that only two coupling constants are non-zero. We thus determine bounds on 
either the 1^ coupling constants of the meson decay (section ||), 



or the coupling constants of the I/-DIS (section 

PaiKji^^ks), Pa{^3lkK2k)^ P'ri^SlkKsk) ^ P(^{^'mk'^'22k) ^ P'T{^32kK3k) ^ Vexp ■ (9) 

There are two special cases where a combination of two coupling constants contribute 
to both lepton flavour violating meson decays and Z/-DIS. Since this is experimentally not 
distinguishable we then consider the combined effects, which also lead to an enhanced 
sensitivity (section |^) 

Fn ■ PA^SlkKlk) + P^ri^SlkKlk) ^ (l - -^tt) ■ PK{>^32kKlk) + P<Ti^32kKlk) < Vexp ■ (10) 

We do not discuss bounds on neutrino oscillation. 

We first continue with a discussion of the NOMAD and CHORUS experiments. We 
then discuss the case of an anomalous meson decay in section ^, the Z/-DIS is treated in 
section ^, and the case with meson decay and u-DlS combined in section ^ In section ^ 
we give a brief summary. Note that throughout this article we present the bounds on the 
^p coupling constants in the mass eigenstate base, in order to avoid model dependent 
results, see ref. [|1^, For theoretical models which predict a hierarchy of ^p coupling 



constants, which we implicitly make use of, see for example |T^, |T 

3 NOMAD and CHORUS 

The NOMAD and CHORUS experiments were designed to search for the oscillation 
of a muon-neutrino into a tau-neutrino through the detection of the charged current 
interaction 

Ur + N ^T + X, (11) 

followed by the decay of the tau in the detector. The search was performed in the muon- 
neutrino beam of the CERN SPS, which provides an intense beam with an average 
neutrino energy of 27 GeV. With a fiducial mass in the ton range (770 kg for CHORUS, 
2.7 Tons for NOMAD), and a number of protons on target of ~ 5 x 10^^, over a data- 
taking period of four years, each experiment has collected on the order of 10® charged 
current muon-neutrino events (i^^*"). Once analysis efficiencies have been taken into 
account, this yields a sensitivity to the presence of tau- neutrinos in the muon-neutrino 
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beam at the level of about one part in 10^ per experiment. The two experiments have 
taken complementary strategies to the detection of tau decay products. 

The iVeutrino Oscillation iVfognetic ZJetector (NOMAD), described in detail in ref. 0, 
is based on an active target consisting of drift chambers immersed in a 0.4 T magnetic 
field. The target is followed by a set of specialized detectors designed to achieve an 
efficient identification of electrons and muons. This design allows to fully reconstruct 
the kinematics of the events, which is used to discriminate between tau decays occurring 
in the detector, and misidentified u^'" and u^'" interactions. A set of detailed analyzes 
which fully exploit the discrimination power of the event kinematics separately for each 
decay channel of the tau have enabled them to put a stringent limit on the muon-neutrino 



to tau-neutrino oscillation probability. The final NOMAD result at 90% CL is ||T7| 

^jrSf^ < 1-63 X 10-^ (12) 

The CERN Hjhnd Oscillation iiesearch Apparatus (CHORUS), described in detail 
in ref. uses a hybrid setup, with a target based on nuclear emulsions, followed by a 
spectrometer system and a calorimeter. This allows the measurement of the momentum 
of the particles emerging from the target region, and the identification of the muons. 

The information of the electronic detectors is used to apply loose kinematic selections, 
which define the data sets on which to perform the tau search in the emulsion stack. The 
decay modes of the tau into a muon or a single hadron are identified by the detection 
of a kink, i.e. a track from the interaction vertex showing a change in direction after a 
short path, of the order of 'ycT^- = (9(1 mm). The latest CHORUS results, based on the 
analysis of a fraction of the events, see ref.|TB|, give no candidate observed, yielding a 
limit: 

P^I^S"^' < 3.4 X 10-^ (13) 

It should be noted that the CHORUS oscillation limit is obtained with a different statis- 
tical treatment than the NOMAD one. If the same treatment were applied, the CHORUS 
result would be roughly PoISi^^ < 2-1 x 10"^, see ref.[||. 



No official combined limit from the two experiments is available. An exercise of 
combining the results, based on the statistical techniques used in NOMAD |T^ , yields 

^oTc^S" < 5 X 10-^ (14) 

We thus use r^exp = 5 x 10~^. 

The two experiments have also obtained limits on the electron-neutrinos oscillating 
to tau- neutrinos, based on the ~ 1% contamination of electron- neutrinos, in the muon- 
neutrino beam. The limits are about two orders of magnitude worse than the previously 
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discussed case, which is why we do not treat electron-neutrinos in this paper. The 
extension is however straightforward. 



n+ or K+) -^fi^ + u 



T 



Starting from eq.(|3), the results in section 2.1 of ref. pO[ lead to 



T^SM -1^11231 ^ -E^' U^J 



with 

2^2 Gf \Vii\ m^(m„ + m^) m 



-^1123 = ^ /- ^ — T— ; 2 ■ (1") 



Gp is the Fermi constant and Vji is an element of the CKM-matrix, m symbolizes a 
mass. Here we have neglected the correction factors due to higher order electroweak 
leading logarithms, short distance QCD corrections, and structure dependent effects. The 
analogous constant to eg. (|T6|) for kaon decay, is obtained from -ft'1123 by replacing 

niK, rud and \'^^^ \^k2 ^ Ki2^^k2- 

With Gf = (0.116639 ± 0.000001) x (100 GeV)-^ = (105.6583568 ± 5.2 x 10"^) 
MeV, = (139.57018 ± 0.00035) MeV, m„ + = (8.5 ± 3.5) MeV, \Vii\ = 0.9750 ± 
0.0008, niK = (493.677 ± 0.016) MeV, m, = (122.5 ± 47.5) MeV, = (21 ± 4) rud, and 



Vul = 0.222 ± 0.004 |]2l[, one gets, using central values. 



|Alr„A3d < ,.,uM(^)\ (17) 



Afci2A3fc2| < 0.0039 



1-F^\ 100 GeV , 



The value of cannot be evaluated with the published numbers. It is in fact a number 
which depends on the neutrino energy, and varies between ~ for the lowest accessible 
neutrino energies to ~ 1 for the highest neutrino energy, as can be seen in Fig.|T|. This 
factor in turn needs to be convoluted with the efficiency of the analysis cuts at a given 
energy, which is also a function of the specific analysis. 

The dependence on the neutrino energy of the global analysis efficiency for CHORUS 
is given in |T^. By convoluting the analysis efficiency with the different components of 



the neutrino spectra shown in Fig. the approximate value for obtained in this way 
is ~ 0.92. We expect that the effective from the combination of the different analyses 
from NOMAD to be lower than this value. In fact CHORUS, for an equivalent limit 
on the oscillation probability, is able to exclude a significantly lower value of Am^, the 
difference of the mass squared of the two neutrino flavours [|l^, |18| . This means [^] that 
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100 150 200 250 300 

Neutrino energy (GeV) 



Figure 1: Spectrum of the muon-neutrino beam at the NOMAD detector, from / |^/ . We 
show the muon neutrino flux originating from different particle decays. The two compo- 
nents of the neutrino spectrum, the one coming from pion and the other one coming from 
kaon decay, are separately shown on the graph. 



the average energy of the neutrinos surviving the analysis cuts is higher for NOMAD than 
for CHORUS. As can be seen from Fig. a harder energy spectrum enhances the fraction 
of events for which the parent meson is a K. We extract a conservative hmit from the 
available information by assuming that the NOMAD and CHORUS results contribute 
with a similar weight to the combined limit. For NOMAD we take F.^ varying between 
~0.1 and ~0.9. The final range for the effective value of in our combined analysis is 
therefore: 

0.5 < < 0.9. (19) 

To be as conservative as possible we use the lower bound when dealing with the pion and 
the upper bound when dealing with the kaons. With eq.(|I^) we therefore obtain 

/ m — ' \ ^ 
|A^*nA322| < 1.4x 10-^ ' I 



< 



X 10" 



100 GeV , 



100 GeV , 



(20) 



Using the product of the bounds on single coupling constants one gets a stricter 
bound on |A']*]^^A3i2| than we would obtain here. 
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In order to determine a bound from Z/-DIS we must compute the cross section 



a 



all Ey 



(21) 



for both the SM process (/ = 2) and the ftp process (/ = 3). g{E''^^) is the incoming 
neutrino energy spectrum (in the lab frame) and is given in ref. [^]. s = {p^^ + PnY = 
2El^^mi\f _j_ jg ^YiQ center-of-mass energy. We perform the above integral numerically. 
The cross section in eq.(|2T|) implicitly contains an integral over the quark distribution 
functions, which we discuss below. Explicitly in terms of the cross sections we have from 
eq.(l) 



rplab 

X E / 9{Elf) ix^u. ^fUN,.-.r+xJs] + X^.c^c aj^^ ,c^r+xjs]) dE, 

rrmi 



'lab 
V 



< 



Vexp- (22) 



E . J 9{Elf) J^] + ,o^,^xJs]) dElf 

rnin. 

We calculate the cross sections for both the SM and the case by first computing the 
partonic cross sections cr^^+q„^fj+g^- With s = {p^^ +PqS^ and t = {p^^ - Paf the 
averaged matrix elements squared in the SM are given by 

= 16G^ |V^,.r (5-m^O (s-mj-m^,), (23) 
(l-^.^+„.cSM^^,.cH = 16 Gl |V,.|2 {s + t- ml, -ml)is + t- ml) . (24) 
In the ^p case, we obtain 



+ Re 



1 

+ 2 



m- 



{s-mli) {s-ml-ml,) 



E 



{\M 



|2\ ^ 



E 

k 
1 

2 



(s + t — — m^i) 



{m\. + m\, - t) {ml - t) 



(25) 



E 



mi 



+ Re 



mi 



+ - 



^ ^fcji ^2fc3 



[s + t — — m^j) (s + t — m%) 

mdr mr {s- ml,) 
:ml+ml-t){ml-t). (26) 



9 



(...) denotes the average/sum over initial/final state spins. The differential cross section 
in the partonic CM-frame is given by 

out 

dcosdC'^ 32 vr s \p,J^^^ ' ^ ' 

where •t}'"^^ denotes the angle between the incoming neutrino and the outgoing charged 
lepton in the partonic CM-frame. In order to obtain the nuclear cross section we have to 
fold the partonic cross section with the pdf's q^{x,t) 

Xmax 1 



da 



•^nmi 1 



d cos 'd'-^^'^ 



{s,t) qZ{x,t) dcos^^^ dx 



{2i 



We use the pdf's given in ref. I^Sj. The Mandelstam variable t is given by0 



(29) 



and y... = y^^/ + ij^l^ — s)^ — 2 m^^fim^^ + s). Furthermore x is the Bjorken scaling 
variable, defined as = xp^. By definition Xmax = 1- 2; is given explicitly in terms of 
the momenta by 



s — m 



X 



N 



2 



A 



'^{Pli + Pqout) "^iV\2 



s — ni 



N 



) -1 



(30) 



We obtain the minimum allowed value, Xmini when the outgoing fermions are at rest in 
the partonic CM-frame 



s — m 



N 



2 m% 



\ 



1 + 



'2{mif + rriq^) tun 



1 



s — m 



N 



(31) 



For each matrix element squared in eqs.(p3D-(p^), we can compute the total cross section 
in eq.(^) using the differential cross section in eq.(P7D as well as the formulae of eqs.(pp|) 
through (0). The SM cross sections as well as the relevant cross section are then 
inserted in eq. (^21) in order to obtain a bound on the ^ coupling constants. 

As stated above, we assume that in any given case only two ^p-couplings are non-zero. 
There are then no interference terms in the matrix elements squared given in eqs.(p5|) 
and ( p6D and the product of the couplings factor out of all the integrals. We thus directly 
obtain a bound on the product of the coupling constants considered.^ 

We determine explicit bounds with the following numerical values rric = (1250 ± 100) 

MeV, TUb = (4200±200) MeV, = (1776.99±0.29) MeV, rup = (938.271998±0.000038) 

^For the sake of generality we shall at the moment not work with = 0, although it is required 
by the parton model. 

•^From now we work with niq^ = 0, as required by the parton model. 
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MeV, m„ = (939.565330 ±0.000038) MeV, |Vi3| = 0.0035 ± 0.0015, \V2i\ = 0.222±0.003, 
\V22\ = 0.941 ± 0.008, \V23\ = 0.040 ± 0.03, We then obtain 



- v/0.43 Xdu + 0.036 Xuc,c V F ilOOGevJ " ^^^^ 

The coefficients of the Xqinqout depend on the pdf's and we determined them numerically. 
Thus each process listed in eq.(^ will have different coefficients. The other processes, 
however, do not lead to improved bounds beyond those published in the literature ||^, 
as we discuss below. For completeness we list the corresponding set of equations in the 
appendix. 

The factor x is 1/0.48 for the NOMAD DIS analyses and 1/0.53 for CHORUS (see 



m- 



k 



ref . and ref . ) . We shall hence work with 

X = 2. (33) 

The standard model z/^*" cross section is dominated by u-rd —>■ t^u. One can therefore 
make the simplifying assumption 

Xdu-- I. (34) 

There is no reliable estimate for the other quark combinations, XqinQout- f^^t the dif- 
ference of the momentum distribution of the quarks in the nucleon, with respect to the 
valence d quark, and/or the presence of heavy (c or h) quarks in the final state signifi- 
cantly alters the kinematic distributions of the final state products. It would therefore 
be necessary to perform a detailed simulation of the experimental analysis to extract the 
efficiency values. Therefore, we only explicitly calculate the case (with Xu^dc = 0) 

|A?„A.„I< 0.0076 (^^5^). (35) 

This is to be compared with the product of single coupling bounds summarized in ref. ||^. 
For k = 1 the product of single bounds is significantly stricter, = 2 is not allowed due 
to the anti-symmetry in the LLE^ operators. For k = 3 the product of single coupling 
bounds is virtually identical to the above bound. However, there the bound on X'^-^^^ 
depends on the squark mass m~ and the bound on A233 depends on the right handed 
slepton mass. If m~ ^ m ~ , our bound will be the best bound. 

The next most promising case for an additional new limit is A^*2iA2fc3, with an in-going 
valence d-quark and an out-going charm quark. Although we cannot quantify Xdc, we can 
evaluate the limit under optimistic assumptions. For Xdc = 1 and XcC'dC = we obtain 



2 



|A-,A2,3|~ 0.006 ^^^) . (36) 
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This is comparable to the product of the single couplings limits for k = 1, which is 
0.002 and somewhat better for k = 3: 0.035. The bounds on the other products are 
not competitive. In the appendix we have summarized the formulae for the bounds from 
Z/-DIS, explicitly retaining Xq„qout,Vexp,X- 

6 (tt^ or K^) — ^ jjP+Ur combined with v^+N — ^ r+X 

Here we use eq. ([IO|) , i.e. we combine DIS and meson decays. For the DIS analysis we 
use eg. (|3^) . For the meson decays we again make use of the results stated in ref . [pOf 



section 2.1. The conservative approach we follow is to set to the smallest value in 
the established range for each of the expressions stated below, and furthermore to set 
XuCdC = Xdc = Xc^dP = 0. We obtain the limits 



|w* w I < 0.117 v/32^ / 

' '''' - v'^+x(0.49 + 0.012 xuodc) llOOGeV, 



/ m~k \ 2 



„ >, , ^ 0.117 ^321^;; / '"S' 



+ x(0.43 Xdc + 0.0015 XcOdo) V 100 GeV; 

^ — I- (38) 

100 GeV J ^ ' 

Note that in the second case we combined jjp + 1/^ with u^ + d t + cj 1^^ + c'^ 

T+dP . Both bounds are improvements compared to products of bounds on single coupling 
constants: respectively 0.0064 and 0.031 

7 Summary and Conclusion 

We have investigated the bounds that can be obtained on the product of -R-parity vio- 
lating coupling constants by employing the neutrino oscillation analyses of the NOMAD 
and CHORUS experiments. They are sensitive to lepton number violating meson decays 
leading to tau neutrinos and to neutrino deep inelastic scattering in which tau leptons 
are produced in the final state. We obtain the following bounds which are more sensitive 
than the product of existing single coupling bounds, and thus present the best bounds 
on these coupling constant combinations 

/ m~2 \2 

lAo*!! A322I < 1.4x10-^ — , (39) 

I 211 3221 _ ^^QQ Q^^j , K J 
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|A^*12^3fc2| < 8.8x10 



-5 / tL 



100 GeV , 



m~2 



/ m~k ^ 

l^m Aoi J < 2.7 X 10"=^ ^ — 

I 32fc 21k\ - I ^QQ Q^y 



(40) 
(41) 
(42) 



For completeness we list the formulae for all processes corresponding to eq.(0) in the 
appendix. Even for optimistic estimates of the efficiencies, the resulting bounds are 
weaker than existing bounds or only a marginal improvement. 
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Appendix 



From Z/-DIS we obtain the following bounds, explicitly retaining XqinQautiVexpiX'- 

\ 2 



^31k ^23fc 



^32A: ^22fe 



'^32k ^23A; 



Afel2 A2fc3 



Afcia A2fc3 



< 



< 



< 



< 



< 



< 



< 



1 



exp 



V0.055 Xsu + 0.028 Xucsc 
1 

v/0.00052 Xbu + 0.0024 XuC^ 



X V GeV 



I 'lexp I 



X V 100 GeV ^ 



1 I / TTl \ ^ 

-L lf]exp I do'^ \ 



X 



v/0.047 Xsc + 0.0036 XcCsC 
1 

V0.00051 Xbc + 0.00039 x^ 
V0.43 Xdu + 0.036 

1 

V0.023 Xsu + 0.036 XuCsC 
1 

VO.OOOlO Xbu + 0.0020 XuCbC 
13 



X V 100 GeV, 



/ Vexp I du 



X V 100 GeV, 

\ 2 



[z. — / m~ 

''lexp I li^" 



X V 100 GeV, 



if]exp I 



X V 100 GeV, 



/ Vexp ( '"^II'' 

X UOOGeV, 
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